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! SUM.HAEY 

A simplified approximate method of calculating the 
unit thermal conductance along an .airfoil -f. a function 
of distance froiu the leading ed^e, "by ug=j of heat transfer 
data for smooth fjlin^ers and smooth f.ah plft^a , is pre- 
sent2d b He.?t tr^n«?er rat:s experimentally obtained by 
several investigators on morals of airfoil? of sections 
E.A.J. 26, K.A.T. 30 f Clarlc 1* , anji NAUA M-6 are compared 
uith results predicted by the use of this method. Calcu- 
lations of the heat transfer rates to he expected from a 
typical full-scale wins are also given. 



INTRODUCE I OCT 



The design of an effective system for the distribution- 
of heat over win-; surfaces to prevent the formation of ice 
requires a knowlodge of the unit thermal conductances along 
such surfaces. 

In order to utilise existing heat -transfer data to 
oalculate these unit conductances! the following ideal sys- 
tem is defined: The leading ed*e of the airfoil is replaced 
by a right circular cylinder with a radius approximately 
equal to the radius of curvature of the leading edge and the 
upper and the lower surfaces of the airfoil are replaced by 
smooth flat plates. The mechanicm of heat transfer along the 
leading edge of the airfoil then corroaponds to that exists 
ing over smooth right circular cylinders, and the mechanism 
along the remainder of the airfoil is postulated to be equiv- 
alent to that existing over smooth flat plates. 
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SYMBOLS 

A area of airfoil equal to chord times span, ft 3 

Op unit hoat capacity 0 f the fluid at constant pres- 

sure , Btu/lb °3P 

0 L lift coefficient defined in F L •= ^C L pAu bo a (dimensionless) 

d diameter of cylinder, ft 

f average unit thermal convective conductance, for 

length L, between airfoil surface and air, 
Btu/hr ft* op 

f unit thermal convoct lire conductance between airfoil 

x surface and air at any point x, Btu/hr fts °F 

f n unit thermal convective conductance botween airfoil 

^ en _ _. _ i _ M 



CP 



surff.ee and air et any angle <P, Btu/hr ft 3 °P 



f c unit thermal convsctivo conductance between airfoil 

^0 surface and air at stagnation point (^=0°), 

Btu/hr ft 3 °7 

I» L lift force defined by P L = bClPAUqo 3 , 113 

g gravitational forco per unit mass, lb/(lb sec a /ft) 

k thermal conductivity of fluid, Etu/hr ft 3 ( 0 j/ft) 

L length of flat plate equivalent to length nlong sir- 

foil surface measured from point of stagnation, ft 

■ 

n exponent of L in equations (5) and (7) (dimension- 

loss) 

q rate of heat transfer, Btu/hr 

b width of span, ft 

t B temperature of p irf oil surface , °T 

3? arithmetic average absolute temperature of fluid 

and airfoil surface, °£ 
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velocity of fluid defined by equation (12) or from 
pressure distribution, f t/soc 

Uco free-stream velocity, of fluid, ft/ sec .. ... 

x length along airfoil profile measured from. the . 

point of stagnation, ft 

x trans loagth along airfoil profile measured from point of 

stagnation to point of the beginning of transition 
of boundary levyor, ft 

x u length along upper surface of airfoil measured from 

point of stagnation, ft 

x L length along lower surface of airfoil measured from 

point of stagnation, ft 

a angle of attack of airfoil, degrees 

Y weight density of fluid at temperature T, lb/ft 3 

1* absolute viscosity of fluid at temperature T, 

lb sec/ft 3 

■p mass density of fluid at tenperature' T, lb sec 3 /ft 4 

T uniform temperature of fluid far from airfoil surface, 

cd -angle between radius through point on cylinder and 

radius through po'int of stagnation measured at 
axis of cylinder, degrees 

Nu Nusselt number (f c d/k) 

pjr Prandtl /number. (3600^0pg/3c) 

Be Reynolds number (uoo&Y/ii g) . (applicable to cylinders) 

.Box Reynolds- number at' distance x from stagnation point 

^ trans Heynol&s number at point of "beginning of transition 
of- boundarjr layar V x trftilB /n g) 
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DISCUSSIOH OP- METHOD' 
Leading Edge 

The magnitude of the unit thermal conductance at the 
stagnation point of a cylinder can "be expressed (reference 
X, vol II, p. 632) sb 



Nu ■ 1.14 Pr°- 4 Ho 0 - 60 (1) 

where 

Nu = k 

M, Cp g 

Pr m 3600 ^ — 

k 



Uo, d V 
He =* * 



The variation of the unit conductance along the forward 
half of a smooth cylinder may be obtained from tho data of 
Schmidt and Vfenner (reference 2). On the basis of these 
data tho unit conductance at any angle cp, moasurod from 
a radius through the point of stagnation -■- that is, at 
qp = 0° — can be expressed approximately as p. function of the 
unit conductance at tho point of stagnation for 00<q><a:30 0 # 
Inspection of these data (reference 2) reveals that as a 
fair approximation 




This expression nay be utilized to calculate the unit 
conductance along tho cylinder, which represents the leading 
edge of the airfoil in the ideal system under consideration. 



Combining equations (1). and (2) gives 
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f c d 



m.1.14 Pro-* Bao.eo ^ -.j^f) 



(3) 



Thus 



Since-, the term ji" -0 * 1 C p 0 * 4 k 0 - 6 for air can be expressed 
approximately as a power function of the absolute tempera- 
ture T, the equation for f 0 becomes, for 0° < cp < 90° f 

% - O.X M ...» e^)— " (1 - Iftl^) <*> 

whore 

d diameter of cylinder 

u^ free-stream velocity 

Y weight density -at temperature 57 

and 

I arithmetic average of temperatures of free atroam 

and of airfoil surfaco3 



Airfcil Surfaces 

The average unit conductance, with length, along a 
flat plate of length I 1b calculated from the equations 
of Oolburn (reference 3, which are usad in reference 4). 

jfor the laminar regime, 



f c a/nL /^m 1» Vvo. bo , v 

Pr a/3 = 0.66(~ j (5) 



3600 u m Y C p • V |* g 



or 



f c = 5.67 (*<>••« n-o- 167 C p °- 333 )"(^Ll) 0,B0 

Again tho properties of air mr*y "bo axprrjssed approximately 
in terms of the absolute temperature T. Then, 

f c = 0.113 T o. so (5L2) 0 -" . (6)- 

3Tor tho turbulent regime , 



o600 



f c rUjj, L Yv-o.ao 
- Pr 2/3 = 0.037 (— ) (7) 



Ins coefficient O.OiS found in references 3 and 4, which 
is tho constant th^t relates the drag coefficient and Reynolds 
number, has beon enacted to agree with that of Goldstein, 
0.. 074/2 (fig, 110, vol. II, pp. 366 and 367 of reference 1).' 



As 



f c = 0.10S /V>.oer c D +0 - 333> ) (^Hl.) 0 -* 0 



then 



f e = 0.655 T°- 396 f^LJL\ 0 - e0 ( 5 ) 
c \^ L o.a5^ 



G?ho unit cor.ductanco .-.t any point x r.lon^ tho flat 
plate iu related to tho ivcrage unit conductance ovar the 
total leiigth L by the equation 
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where f c Is the average conductance for the length x 
from equation (8), and n is the exponent of L in equa- 
tions (5) and (7) (reference H). The unit conductance 
along the flat plate, which in the ideal system is equiv- 
alent to the portion of the profile beyond the leading edge, 
may he expressed "by 



f 0 ° O.O562 T° * B0 
for the laminar regime, and hy 

f c « 0,52^ T 0#S98 
c x 

for the turhulent regime 

where x is the distance alonp the flat plate measured 
from the point of stagnation of the airfoil. The unit 
thermal conductance f c ^ as a function of u m V/x and 

of u m V/x 0 ' 36 for use in equations (10) and (11), re- 
spectively, is shovn in figure 1. 

Because of the circulation around the actual airfoil, 
the velocity of the fluid is greater along the upper sur- 
face and smaller along the lover surface than the free- 
stream velocity u^. In reference 5 it is demonstrated 
that, as a first approximation, th* velocity to "be used 
in equations (10) and (ll) for the low^r and upper sur- 
faces of the airfoil is 




^ (l * ) (12) 

\ 4 cos a/ 

in which C L is the lift coefficient of the airfoj.1 de- 
fined hy P L = ^ C L pAuoo a and a is the angle of attack. 

The positive sign in equation (12) is used in the compu- 
tation of the velocity along the upper surface and the 
negative sign is used for the velocity along the lo^er" 
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surface* The velocity distribution derived from pressure 
. measurements along the airfoil may also he used for the 
evaluation of u„ in equations (10) and (11). This pro- 
cedure yields slightly higher values of f c near the 

leading edge on the upper surface tut the average unit 
conductance, with length, does not differ greatly from 
that obtained from equation (12). The value of f Q ^ on 

the lower surface, where the flow is usually laminar over 
a large part of the airfoil, is not so greatly influenced 
"by the magnitude of u m which entors into equation (10) 
as the 0.50 power. 



Transition "between Laminar and Turbulent Boundary Layer 

In order to calculate the unit conductance of an 
actual airfoil at a givau angle of attack a, an equiva- 
lent cylindor is substituted for the leading edge and the 
local unit conductance is computed up to an angle of ±90° 
by use of equation (4) and the froe-strcam velocity 
The variation of tho position of the stagnation point with 
the angle of attack is neglected in this calculation be- 
cause its change of position is small compared to the chord 
of the airfoil. 

The unit conductance for the portion of the airfoil 
following the equivalent cylinder is calculated from the 
f let-plate equations (10) and (11), by using the distance 
along the profile measured from th? point of stagnation 
x and the corrected velocity u m . The value of the abso- 
luto temperature T in equations (4), (10), ?.nd (11) is 
dof ined as the nri thmc t ic mean of the fluid and the air- 
f oil temperatures • 

An examination of equetlon (10) indicates that, for 
tht) region of laminar flow, the m%gr.itudo of f actually 
is indopondent of the* arithmetic morn absolute temperature 
T. This result is due to the fact thai; the proportics of 
air k°- ae7 jj.-o.i67 Cp o.333 aro proportional to the 0.50 

power of T while the woight density y f which is inversely 
proportional to the absolute temperature T, enters with the 
0.5C powor in equatioa (6). The magnitude of f<, in equa- 
tion (4) is practically indopondont of T, since- tho prop- 
erties of air - that is, k, n, and C - are a function of 
the 0.49 power of T. Jor tho rogion^of turbulent flow 
(equation (11)), a similar procedure revoals thr»t the local 
unit convective conductance is proportional to tho -C. 60 
powor of T. 
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An approximate picture of the flow conditions in the 
"boundary layer along a flat plate is shown in figure 5 of 
reference 6. 

The poin.t of - tran si t ion of the boundary layer from 
laminar flow to turbulence is indeterminate* The position 
of this point is a function of the free-stream velocity and 
the corresponding turbulence , the angle of attack , the shape 
of the airfoil i the surface roughness, and other variables* 
The* point of transition along the upper surface may be esti- 
mated roughly, however, by setting . B e trans a "^^^^^ ~ 

(reference 1, vol I, p. 326) equal to a value between 

5 x 10 and 5 x 10 a . These magnitudes of Be tranB refer 

to the point at which turbulence begins. The rate at which 
turbulent flow becomes fully developed may bo estimated 
from the behavior of the experimental values of fc z « Ehe 
position of the point of transition along the lower surface 
is also unknown, although some evidence indicates that the 
boundary layer along tho lower surface of models is laminar 
for a considerable part of the chord at angles of attack 
greater than 0°, 



Rate of Heat Transfer from Airfoil 

Figure 2 illustrates th* distribution of the point 
unit thermal conductance f Cx obtained by the irethod pre- 
viously outlined, Tho lower surface would hpve n distri- 
bution similar to that shown along tho upper surface but 
would yield somewhat lowor magnitudes of f Cx because of 

lower values of , and the transition point would occur 
at a greater distance from the leading edge. 

The maximum average unit conductance over the whole 
airfoil occurs when the point of transition Is close to 
the nose. As a conservative estimate, tho heating system 
s.hould bo designed for this position of the transition 
point, since any disturbance, such as initial ice formation, 
would tend to move the transition point forward. 

The r^.te of he?.t transfer from the whole airfoil may 
bo calculated by the equation 




x 



/ 
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where 



span 

heated length, measured along profile on upper 
surface 

heated length, measured along profile on lover 
surface 

difference in temperature between airfoil sur- 
fnoe and ambient fluid at point x measured 
from leading edga- 



DISCUSSION- OF EXFESIMMTAL DATA 
E.A.J. 26 Model 



Observations of the thermal characteristics of an 
E.A.F. 26 airfoil (reference 7) are given in figure 3, 
which also shows the pradictod values of f~ as a funo- 

^X 

tion of the distance x along the profile. The tests 
were performed on a 6-inch airfoil model, which was heated 
electrically by means cf platinum resistance strips placed 
spanvise on both upper and lower surfaces. 

The predicted unit conductances for the nose were set 
equal, as proposed horeiuj to the unit coy.ductn.nces along 
a cylinder having a radius equivalent to thti radius of 
curvature cf the lsadlng odge. The variation of the unit 
conductance with distance along the upper p.nd lower surfaces 
was then obtained from equations (10) and (11) by using the 
corrected velocities u m . 

The experimental result for the unit conductance at 
the ncse is always lower than the predicted vriue at the 
point of stagnation because the measured vcJuo is an average 
over a heated strip which includes jnore- surface than the 
point of stagnation. Th3 position of this point is a func- 
tion of the angle of attack, but the predicted values are 
based on a point of stagnation fixed p-t the leading edge. 

On the lower surface of the airfoil at a = -0.9° 
(fig. 3), the experimental results can bo predicted for 
a given position of -the transition point. The experimental 
data indicate that this point moves toward the leading edge 



(*.-T.) x 
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as the free-stream velocity is increased. At a = 3- ?0 
(fig. U) f the flov is apparently laminar along most of 
the length of the surface and the unit conductance can 
"be closejly predicted... The agr_eerient of -predicted and 
experimental results "becomes "better as the velocity 
is decreased and as- the angle of attack a is increased, 
"because the laminar layer is more stable at small magni- 
tudes of u m . Owing to the circulation around the air- 
foil, the local velocity u m , calculated from equation 
(12), is less than the free-stream velocity tx^ and de- 
creases as a increases. The -predicted results for the 
lower surface are inappreciably affected "by the use of 
the tabulated values of the fluid velocity just outside 
the boundary layer. 

Along the uoper surface also the unit conductance can 
be predicted if the transition point is known. As might be 
expected, the experimental results for or. = -0.9° show that 
the transition point occurs far back along the trailing edge 
of the 6-inch model and that it moves forward as u^, increases, 
For a = 3.2° and a = 7° (figs. ^ and 5), the calculations 
for the ideal system do not predict th*» v«*ry high values of 
the- unit conductances that' apparently occur in the region 
of transition from laminar to turbulent flnv. This deviation 
between the predicted and the experimental results could be 
made smaller if the distance along the equivalent flat x>late 
x were measured, not from the point of stagnation, but from 
some point farther back along the plate. The deviation may 
be ascribed to the great curvature of the airfoil in the re- 
gion of the leading edge, ^hich may invalidate substitution 
of a flet plate for the airfoil in that region, and the use 
of u ffl . The change in x would not appreciably aff9ct the 

predicted results along the rear half of the chord but would 
have a great effect upon the slope within the first 70 per- 
cent of the chord, since f_ varies as x""°' a . The slopes 

c x 

of the theoretical curves of f c ^ against x approach small 

and almost constant magnitudes with Increasing values of x. 
The experimental values along the upper surface rapidly dimin- 
ish at the trailing edge; the experimental values along the 
lower surface rapidly' rise and, at th« end of the model air- 
foil, approach the same value as those along the upper surface. 
The previously mentioned deviations between the predicted and 
the experimental results near the region of transition can 
also be mp,de smaller (fig, **(c)) if the tabulated velocity 
distribution is used in equations (10) and (ll) instead of 
the value of u^ derived from enuation (12).. This method. 
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however, does not predict the very high value of the unit 
conductance which occurs at the point where the "boundary 
layer is apparently fully turbulent. The corresponding 
result near the trailing edge is the same an that vhen 
the value of u m from equation (12) is used. 

R.A.F. 30 Model 

The experimental results on a 2^-inch model of an 
B.A.J 1 , 30 airfoil at a = 5. 3° are compared with predicted 
values in figure 6. The heated surface, located at the 
nose, comprised RT-proxim/it ely 25 percent of the chord for 
the upper surface and 1^ percent of the chord for the lower 
surface. The total heat loss of the entire heated surface 
was determined in the experimental work (reference g) for 
several angles of attack at velocities ranging from Uo feet 
per second to 130 feet per second. 

The predicted variation of the unit conductance was 
determined according to the methods proposed herein, after 
arbitrarily locating the point of transition at which tur- 
bulence begins. Tor the up-rer surface this -point wes taken 
at different velocities for R fi trans from ?0»000 to 300,000, 
which corresponded in each case to transition at a point 
close to the nose. For the lover surface these points of 
transition were taken at He trans from 30,0^0 to ^00,000. 
These predicted valueq of the unit conductances wore averaged 
with respect to distance along th° airfoil over the entire 
length of the heated section (>oth upper and lower surfaces) 
fnr each velocity. Th6 averaged predicted values nre about 
20 percent lower at all velociti°s vhon u m *-as determined 
from equation (12) and about 10 percent lover vhen th° ve- 
locity distribution dorived from Treasure measurements along 
the airfoil was used for the values of u m than the corre- 
sponding experimental results. 

■ Clark T Model 

The data for a Clark Y airfoil '"ere taken on a 10-inch 
model (reference ?), vhich wps heated in four separate sec- 
tions along the chord, including th5 leading and the trail- 
ing fldpes. The predicted vp.lues were determined for one 
velocity and for two different angles of attack. Data and 
predicted values are compared in figure 7» The predicted 
unit conductance over the upper and lower surfaces included 
points of transition for which th r values of the Reynolds 
number were. 2.3 x 10 4 at a = 6° and 3.2 x 10 4 it a = 0° 
for the up^er surface and 1.5 x 10 B at r = 6° and 1.3 x 10 B 
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nt a » 0 for the lower surface. The predicted values 
again appeared to he, at the moot, ahout 20 percent lower 
than the experimental results averaged with respect to 
length for each section. The values of u m were calcu- 
lated from equation (12). 



NACA K-6 Model 

Experimental data taken on a 10-inch model of an 
HAG A M*-6 airfoil (reference 10) are presented in figure 
8 for . a = 4° and for free-stream velocities ranging from 
96 feet per second to 303 feet per second. The model was 
heated at the nose and at three contiguous sections on the 
upper surfaco and on the lower surface. The heat transfer 
from each heated section was determined and is shown in 
figure 8 as constant values over each section. The expert 
imentally determined conductances are superposed upon the 
predicted distribution of conductance along the airfoil. 
If the flow over tho entire model were considered to ho 
laminar along most of the lover surface, tho experimental 
results could bo predictod closely. The measured value* 
of the conductance along the upper surface clearl7 indicate 
the increase in tho value of f c when the "boundary layer 
"becomes turbulent, values for t£o unit conductance for tho 
case vhon the no&e section was &ot heated arc also shown 
in figures 8(h) and 8(d), Th; fact that the heat transfer 
from the other soctiona is higher when the noso is not haated 
may "be ozplained by a delayed transition "of tho boundary 
layor from laminar to turbulont flo*r due to the removal of 
the inital disturbance; thnt ie f the heated noso soction. 

It should be kept clearly in mind that, for a model, 
tho distribution of the conductance depends a great deal 
on the position of the transition point; but, in the case 
of a full-sise airfoil, the transition region is a small 
part of tho total area usually hoated and for this -reason 
its position need not be accurately known. 



SrTJMEEICAL EXAMPLE 



As a further illustration of the application of the 
proposed mothojl for computation of the distribution of the 
unit conductance along a full-size airfoil, an oxamplo is 
presented. 



The distribution of the unit conductance along an 
actual airfoil with a 6-foot chord find a profile which 
approximates that of the UAOA 23012 is shown in figure 9 
for a s Qo and for a free-stream velocity equal to hko 
feet per second (300 mph ) . The transition point for the 
"beginning of turbulence for these conductances was taken 
at a value of He trana of approximately 0, which is at 
the stagnation point, and at 500,000. 

If the airfoil were heated for a distance of 2 feot 

along "both the upper and the lower surfaces, the predicted 

average unit conductance for this section, when values of 

u n calculated fror equation (12) are. used, would "be only 

10 percent lower at a free-streain velocity of 300 miles 

per hour and 23 percent lot'er at a fre*-strenm velocity of 

150 ni>a per hour (see fi£. 9), when Ee trans is taken 

to he 500,000 than the unit conductance when Re, is 

trans 

at the leading edge - that is, whrn R^trans =0, 

These calculations demonstrate the fact that for a 
full-site airfoil the exact position of the transition 
point "becomes loss important at higher velocities *nd for 
larger heating sections, since a gre*t variation in tho 
position of this point was arbitrarily chosen - that is, 
He trans at 0 and at 500,000. 

The distribution of the unit ronfluctancp when the 
value of Urr is taken from pressure rceasur pment s along the 
airfoil is nlso shown in fig-ure 9- The average value of 
f c , vith length along thp airfoil, is ahout 7 percent 
higher than that if u r is evaluated fron equation (12). 

Figure 10 reveals the predicted values of the h^at 
loss from an KACA 23012 airfoil for a difference of 70° F 
"between the temperatures of the fluid and of the surface 
at velocities ranging from 150 to 300 miles per hour. The 
magnitude of the -weight density Y ras evaluated at atmos- 
pheric pressure and at 0° F, which is the arithmetic average 
of the temperatures of the free str eam:-.*nd of the airfoil 
surface. The values of u w*ro calculated from equation 
(12). 

For the case of an actual tapared wing, the value q/h 
should "be computed for each foot of span "because of the 
variation in profile and chord length. The total rate of 
heat transfer would then "be the sum of these q/h contri- 
"but ions • 



15 



OOKCLUSIOHS 



1« The variation of the unit thermal conductance 
along a wing can be estimated fairly accurately by an 
approximate method from known heat-transfer data on 
smooth cylinders and flat plates. 

2m An exact computation of the distribution of the 
unit conductance along a wing requires the determination 
of the- position of the transition of the boundary layer 
from laminar to turbulent flow. Such knowledge is nec- 
essary for the proper design of the heat— .distribut ion 
system in the leading edge. 

3. A precise agreement between the predicted results 
and experimental data taken on small airfoil models is not 
to be expected because the transition pdint may Occur over 
a large percentage of the chord in the case of models, 
depending on the conditions of flow; but this change in 
position comprises a email part of the chord in the case 
of full-size wings. 

4. The method proposad herein is satisfactory, how- 
ever, for a conservative estimate of the thermal capacity 

of a de-icing system. The Accuracy of the method ia greater 
for high velocities and for heating sections that extend 
ovor a larger part of the profile, She velocity along tho 
airfoil can be obtained from s tat ic-prosaur o measurements 
or, if th^se dntft are not available, can be estimated from 
an equation givon herein. 



University of California, 
Berkeley, Calif. 
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FI9.IO- Predicted rate of heat transfer as a function of distance along the airfoil surface 
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Figs. 3a, b 
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Fig. 3o Unit thermal conductance as a function of the distance along the 
a irf oil surface. 
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Fig. 3b Unit thermal conductance as a function of the distance along the 
airfoil surface. 
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Figs. 3c , 4a 
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Fig. 3C Unit thermal conductance os a function of the distance along 

♦ he a i r f o i I s u rf a ce . 
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Fig. 4CL Unit thermal conductance as a function of the distance along the 
airfoil surface. 
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Figs. 4b , c 
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Fig. 4- b Unit thermal conductance as a function of distance along the 
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Unit thermal conductance as a function of the distance along the 
airfoil surface. 
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Figs. 5a , b 
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FIG. 5b UNIT THERMAL CONDUCTANCE AS A FUNCTION OF THE DISTANCE 
ALONG THE AIRFOIL SURFACE. 
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Fig. &a Unit thermal conductance as a function of the distance along the 
airfoil surface. 
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Fig. 6 Predicted and miaiurid oviragi unit thermal 
conductance of RAF 30 Airfoil as a function 
of velocity. 
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Fie* 7 Unit thermal conductance as o function of the distance along 
the airfoil surface. 
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Figs. 8b ;c 
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Fig. 8b Unit thermal conductance as a function of the distance along the 
airfoil surf a ce. 
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Fig. Sc Unit thermal conductance as a function of the distance along the 
a i r f o U surface. 
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Figs. 8d, 9 
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Fig. 8d Unit thermal conductance as a function of the distance along the 
airfoil surface. 
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Fig. 9.- Predicted unit thermal conductance as a function of the distance 
along the airfoil surface. 
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